Effects of Electronic Correlations on the Thermoelectric Power of the Cuprates 
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We show that important anomalous features of the normal-state thermoelectric power 5* of high-Tc materials 
can be understood as being caused by doping dependent short-range antiferromagnetic correlations. The theory 
is based on the fluctuation-exchange approximation applied to Hubbard model in the framework of the Kubo 
formalism. Firstly, the characteristic maximum of S as function of temperature can be explained by the anoma- 
lous momentum dependence of the single-particle scattering rate. Secondly, we discuss the role of the actual 
Fermi surface shape for the occurrence of a sign change of 5 as a function of temperature and doping. 



PACS numbers: 71.27.+a, 74.25.Fy, 74.20.-z, 74.72.-h 

The thermoelectric power (TEP), S, of the high-Tc materi- 
als in the normal state exhibits anomalous features that 
are not well understood at present. For example, its temper- 
ature dependence shows a characteristic maximum that can 
be found in all optimally and underdoped cuprates, which is 
in contrast to the conventional linear-T behavior of a weakly 
correlated Fermi liquid. Furthermore, a better understanding 
of the doping dependence, which has been shown to be uni- 
versal for a large class of high-Tc materials [|l]], is needed. 

Experimentally, the TEP is positive for the underdoped 
cuprates. It increases rapidly as a function of temperature, 
reaches a maximum S* at a temperature T* and falls off al- 
most linearly with temperature. The size of the maximum S* 
and the value of T* decrease rapidly with increasing doping, 
while the slope for the fall-off at T > T* is almost indepen- 
dent of the doping concentration. The TEP of nearly optimally 
doped materials shows a similar behavior below room temper- 
ature, but the over-all size of the TEP is reduced. Furthermore, 
the TEP for most of the optimally doped samples changes sign 
at approximately room temperature For overdoped sam- 
ples, the behavior of the TEP is very different: It is negative 
and decreases linearly with increasing temperature. At present 
there are only two classes of high-Tc materials that do not 
follow this generic trend, namely La2-xSrxCu04 (LSCO) 
and YBasCuaOr-a (YBCO). For LSCO the TEP remains 
positive in the overdoped regime, where it exhibits a maxi- 
mum of decreasing height for increasing hole concentration 
iQ], whereas YBCO shows a negative TEP in the overdoped 
regime but with a positive slope at temperatures T > T* [|l |2|] . 
Recently Bernhard and Tallon [|[] presented strong experi- 
mental evidence that the chain contribution is responsible for 
this non-generic TEP of YBCO while the contribution of the 
Cu02 planes follows the generic behavior 

Theoretically, important open problems concerning the 
characteristic behavior of the TEP remain, although there have 
been several attempts to explain this behavior, ranging from 
the van Hove scenario to phonon drag effects [^. The 
most important open questions are the physical origin of the 
characteristic temperature scale T* and of the doping depen- 
dence of the sign of S. 



In this paper, we demonstrate that short-range antiferro- 
magnetic correlations are responsible for salient features of 
the temperature and doping dependence of the TEP. In par- 
ticular, we show that for underdoped systems the occurrence 
of a pseudo-gap in the single-particle excitation spectrum and 
for optimally doped systems the pronounced increase of the 
scattering rate near the "hot-spot" regions (near (tt, 0)) are es- 
sential. As was proposed by Stojkovic and Pines [^, the latter 
phenomenon is also of importance for the description of the 
temperature dependence of the Hall coefficient. Furthermore, 
we discuss the influence of the shape of the Fermi surface of 
LSCO and YBCO and of the bilayer coupling. 

We note that antiferromagnetic spin correlations have re- 
cently also been argued to be intimately related to the normal- 
state pseudogap in underdoped cuprates and to a sizeable de- 
formation of the quasi-particle band structure and Fermi sur- 
face as a function of doping [p|jlO|]. These ideas have been cor- 
roborated with evidence from Quantum-Monte-Carlo (QMC) 
simulations in a recent work by Preuss et al. [[ll]]. These latter 
results also provide numerically, in principle exact support for 
the relevance of antiferromagnetic correlations that are found 
to also be crucial for the thermoelectric power in the theory 
presented here. 

The thermoelectric power S is determined by the relation 
between the electrical current j and the temperature gradient 
VT: j = CTdc {'El — S VT), where is the dc-conductivity 
and E an applied electrical field. In order to account 
for the strong electronic correlations and the pronounced 
short-range antiferromagnetic correlations of the cuprates, 
we use the two-dimensional single-band Hubbard model 
with e(k) = —2t [cos kx + cos ky) — At' cos k^ cos ky — 
2t" {cos{2kx) + cos(2fcj,)) for the microscopic calculation of 
S. We set t = 0.25eV and U = 4t. The values of the longer- 
range transfer integrals will be given below. The transport co- 
efficients (Tdc and S are calculated within the Kubo-formalism 
as 



Cdc = e / du) {—dj{uj) / du) a {lo) 
e 



5* = -- 



dw {~df{uj) I duj) uj a{ij:) 
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with the differential conductivity 



(3) 



Eqs. ^ and (||) are obtained following, for example, the 
derivation of Moreno and Coleman [ p^ and approximating 
the full vertex function by the velocity fk.x = ^x^O^)- Here e 
is the electronic charge, N the number of sites, /(cj) the Fermi 
function and yl(k, a;) = — ilmG(k, uj) the spectral function. 
The single-particle propagator G(k, oj) is calculated within 
the fluctuation-exchange approximation (FLEX) |]l3|]. This 
approximation includes the interaction of the electrons with 
charge and spin fluctuations. We employ a real-frequency ap- 
proach to this approximation, which directly yields A(k, oj) 
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FIG. 1. Thermopower 5* for the LSCO system as a function of 
temperature for different doping concentrations. The inset shows S 
of the tight binding model for the same doping values. The triangles 
indicate the high-temperature limit given by the generalized Heike 
formula. 



In Fig. |l|, we show the TEP for the single-layer LSCO sys- 
tem as function of temperature for different doping levels. 
Here LSCO is described by t' — t" — to ensure that the 
Fermi surface is closed around (0,0) [p^. For all fillings, 
we find the interesting result that our theory yields a positive 
TEP which decreases with increasing hole concentration and 
shows pronounced maxima at about 150K for the two lowest 
doping levels. For the higher doping level (20%), the maxi- 
mum is very broad and shifts to higher temperatures. This is 
very similar to the experimental observation in LSCO [| 16 1. 
Note that also the absolute values of the calculated TEP are in 
good agreement with experimental data [Q] . 

In order to understand the role of electronic correlations for 
the T-dependence and sign of S and to relate them to pre- 
vious theories, it is at first of interest to contrast them with 
the TEP of the noninteracting system (U = 0). Here S is ob- 
tained from a tight-binding model calculated in the relaxation- 



time approximation (the U = calculations along the lines 
of eqs. (|l]) to ^ yields a vanishing TEP) with a momentum 
and frequency independent relaxation time, see the inset of 
Fig. |l} In this approach, S is negative and proportional to T 
in the relevant temperature range {T < 600 K). In the high 
temperature limit, the TEP deviates from this linear behav- 
ior for T ^ t and finally approaches the values given by a 
generalized high-temperature Heike's formula [|l7|]. Such a 
U = theory, in which no details of the relaxation or scatter- 
ing rate are included, still fails to describe experiments even 
when a more sophisticated phenomenological band structure 
with a pronounced flattening of the dispersion near the Fermi 
level [ [l8| ] as found in photoemission experiments [ [l^ , ^0| ] is 
included. Thus, the strong qualitative disagreement between 
the above U = results and our FLEX results for the inter- 
acting system implies that for an understanding of the experi- 
mental data it is necessary to take into account the momentum 
and frequency dependence of the scattering rate. In order to 
demonstrate in more detail the role of electronic correlation 
effects for S{T) in our approach, we now show (i) that our 
results for the sign of S are caused by the anomalous k de- 
pendence of the scattering rate 1/T(k) = — Im E(k, = 0) 
(where cu = refers to the Fermi energy) due to antiferromag- 
netic correlations and {ii) that the temperature dependence 
of S is determined by the characteristic temperature scales of 
these correlations. 
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FIG. 2. Differential conductivity (solid lines) and density of 
states (dotted lines) for two doping levels and T — 200 K. All for 

U = At. 

It is instructive to consider the differential conductivity 
cr{uj) of Eq. (H), since the TEP, from Eq. (^, is a measure 
of the asymmetry of a{u!) with respect to cj = for energies 
of the order of ksT. Thus the TEP is very sensitive to excita- 
tions of states which are close to but not directly at the Fermi 
surface. 

InFig.|(a), we show cr(aj) and the density of states g{LLj) = 
^0^1 ^) for ^ doping concentration of 20%: Although 
g{uj) is peaked for > (as in the uncorrelated case), the dif- 
ferential conductivity (j{u!) exhibits its maximum for u; < 0, 
yielding the positive value of S. This behavior can be un- 
derstood if one takes into account that in Eq. (j|) momentum 
states with small velocity Vk and large line width (large scat- 
tering rates) are suppressed. Antiferromagnetic spin fluctua- 
tions lead to large scattering rates for momenta near kp + Q 
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(kp Fermi momentum, Q = (tt, tt)) which, for LSCO, are 
everywhere outside of the Fermi surface (see Fig. ^a)). 
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FIG. 3. Fermi surface, (a,b), and spectral function, (c,d), for 
the LSCO and YBCO parameters, respectively. In (c) and (d) dark 
(white) areas correspond to large (small) spectral weight. 

Thus the scattering rate of unoccupied states (holes) near 
the Fermi surface is larger than the scattering rate of the cor- 
responding occupied states (electrons). This suppresses the 
contribution of the holes (at w > 0) to (j{u!), which as a result 
has a maximum at a; < 0. This effect is illustrated in detail 
in Fig. 3(c), where we show the spectral function A{k, uj) in 
the vicinity of = 0. It can be clearly seen that the spectral 
weight is more sharply peaked for occupied than for unoccu- 
pied states around (tt, 0). The suppression of states near (tt, 0) 
due to the low velocities near the van Hove singularity also 
contributes to the absolute magnitude of S but does not deter- 
mine its sign and order of magnitude. This demonstrates that 
the thermoelectric power is determined by the combined mo- 
mentum and energy dependence of the scattering rate, making 
transport for holes less coherent than for electrons, yielding 
S > 0. The transition to conventional behavior occurs at tem- 
peratures for which the antiferromagnetic correlation length 
becomes smaller than the range of the hopping matrix (deter- 
mined by t, t' and t"), such that no pronounced momentum 
dependencies in the scattering rate occur [pHTTI]. For decreas- 
ing doping concentration, the strength of the antiferromag- 
netic correlations increases, leading to the appearance of a 
pseudogap in the DOS as a precursor of the antiferromagnetic 
phase transition | p^ , pl] |. The FLEX calculation for U — At 
[ pX| ] yields results for the pseudogap qualitatively similar to 
those obtained in numerically exact QMC simulations JTl]]. 
The density of states showing the pseudogap is presented in 
Fig. |^(b). The asymmetry in g^uj) caused by the antiferro- 
magnetic fluctuations is even more pronounced in a{uj) such 
that these coiTelations cause in addition to the positive sign 
a pronounced increase of the magnitude of S. Even more 
interesting, the T-dependence of the antiferromagnetic fluc- 
tuations, which also manifests itself in the pseudogap [|l0|]. 



determines the anomalous T-dependence of the TEP in the 
under- and optimally doped systems of Fig. |l| Note, however, 
the weak doping-dependence of T* in our theory. We believe 
that this shortcoming is related to the difficulties in the theo- 
retical description of the single-particle excitation spectrum of 
underdoped systems rather than to the approximate treatment 
of the transport problem. We believe that a better description 
of the detailed momentum dependence of the pseudogap, in 
agreement with recent angular resolved photoemission exper- 
iments 1 19 1 and QMC results niin is essential here. 
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FIG. 4. TEP for the YBCO parameters as a function of tempera- 
ture for different hole dopings (inset: doping dependence). 

The inset in Fig. ^displays our results for the doping depen- 
dence of the planar contribution to the TEP of YBCO which 
is also generic for other cuprates, using t' ~ —0.38t and 
t" = — 0.06t such that the Fermi surface is closed around 
the (tt, tt) point (see Fig. ^(b)). In contrast to the LSCO sys- 
tem and in agreement with the general trends in the experi- 
ments, the thermopower changes its sign as function of the 
doping concentration. In order to understand the origin of this 
behavior, we plot in Fig. ^ the TEP as a function of temper- 
ature for different doping values. As in LSCO, a maximum 
of the TEP occurs for small doping and small temperatures. 
However, for increasing T, the sign of S becomes negative. 
For larger doping, S decreases until it becomes negative at all 
temperatures. In agreement with the experimental observa- 
tion [^], we find that the negative slope of the TEP for large 
T is almost doping independent. These phenomena can also 
be understood in terms of the pronounced momentum depen- 
dence of the scattering rates caused by antifeiTomagnetic spin 
fluctuations. We find that for low doping, the scattering rates 
are maximal for momentum states near kp + Q. However, 
as can be seen in Fig. |]^b), these states are solely outside the 
Fermi surface for kp + Q near (f, -7), while they are within 
the Fermi surface for kp + Q near (tt, 0), in contrast to the 
LSCO system. If these two regions in momentum space were 
of equal importance for the suppression of near Fermi surface 
states in <j{uj), they would cancel each other, leading to a neg- 
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ative S, as found for f/ = 0. However, the transport properties 
are dominated by the quasi-particles near the diagonal of the 
Brillouin zone |21 1. This can also be observed in the spectral 
function plot for YBCO in Fig. |(d), in which the states near 
{tt/2, 7r/2) are much more coherent than the ones near (tt, 0). 
For kp near {tt/2, 7r/2), the holes near the Fermi surface are 
closer to kp + Q than the corresponding electron states and 
it follows that S* > for sufficiently strong antiferromagnetic 
correlations. Finally, due to the more subtle competition of 
different regions in momentum space for the YBCO param- 
eters, the weakening of antiferromagnetic correlations for in- 
creasing doping or temperature leads to the more rapid (as 
compared to LSCO) transition to negative values of the TEP, 
shown in Fig. Q For a quantitative comparison of our results 
in Fig. with experiment, it should be noted that the mag- 
nitude of S as well as the doping level at which S changes 
its sign for low T are smaller than in experiment. One rea- 
son for this quantitative difference might be the neglect of the 
bilayer coupling in YBCO which was shown to enhance the 
in-plane antiferromagnetic correlations [ |2^ , ^3| ]. By taking bi- 
layer coupling into account via an interlayer hopping t±, we 
solved the FLEX equations and determined the TEP. Note, 
that t± = 0.4t is a reasonable value for the interlayer hop- 
ping in YBCO [p3[]. The results are also shown in the inset of 
Fig. ^ for different values of t± and demonstrate that indeed 
a considerable increase of the absolute magnitude of S occurs 
due to inter-layer spin correlations. 

Concerning the influence of the bandstructure parametriza- 
tion used for YBCO, it might be argued that the velocity Wk.a 
on the Fermi surface depends more strongly on the direction 
of k than for other parametrizations (e.g. t' = — 0.45i and 
t" = [|]]). However this is not essential for our results, since 
the agreement with experiments would be even better if we 
were to neglect this k dependence, which suppresses S. We 
have used this parametrization since it yields better agreement 
of the calculated spin susceptibility with the spin-excitation 
spectrum deduced from NMR-experiments [^J] . 

In summary, we have shown that the anomalous features of 
the in-plane thermopower can be understood as consequences 
of short-range antiferromagnetic correlations. The opening 
of a pseudogap and the anomalous momentum dependence of 
the single-particle scattering rate provide an explanation of the 
general trends of both temperature and doping dependence of 
the in-plane thermopower. This demonstrates that the thermo- 
electric power is a sensitive probe of the anomalous nature of 
the low-energy excitations, caused by the antiferromagnetic 
correlations of the cuprates. 
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